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The two-pion correlation function can be defined as a ratio 
of either the measured momentum distributions or the nor- 
malized momentum space probabilities. We show that the 
first alternative avoids certain ambiguities since then the nor- 
malization of the two-pion correlator contains important in- 
formation on the multiplicity distribution of the event ensem- 
ble which is lost in the second alternative. We illustrate this 
explicitly for specific classes of event ensembles. 

PACS numbers: 25.75.-q, 25.75.Gz, 25.70.Pq. 



I. INTRODUCTION 

Two-particle Bose-Einstein (BE) interferometry (also 
known as Hanbury Brown- Twiss (HBT) intensity inter- 
ferometry) as a method for obtaining information on the 
space-time geometry and dynamics of high energy colli- 
sions has recently received intensive theoretical and ex- 
perimental attention. Detailed investigations revealed 
that high-quality two-particle correlation data constrain 
not only the geometric size of the particle-emitting source 
but also its dynamical state at particle freeze-out 

Two different definitions of two-pion correlation func- 
tion are employed in the literature [p]-p^ . The first starts 
from the measured invariant i-pion inclusive distribution 



which is normalized via 
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to the ith order factorial moment of the pion multiplicity 
distribution, and defines the two-particle correlator as 



Nl{pi)Ni{p2) 



(3) 



The second definition instead employs the normalized i- 
pion production probability 



P,{pi, . . . ,Pi) 
and defines 
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It follows that 

C"{PI,P2) 



{n{n — 1)) 



-yC^(pi,P2). 
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Recently, Miskowiec and Voloshin argued that the 
first definition is preferable because it is based directly 
on measured quantities and it is consistent with the often 
used theoretical expression 



C{q,K) = l 



1/ d^xSix,K) 



Jd^xS{x,Pi) Jd^yS{y,p2) 
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Here K ^ {pi + P2)/2, q ^ pi ~ P2, q" = Ep^ - Ep^, 
and S{x, K) is the emission function of the source. In 
this paper we will stress that the first definition has the 
additional advantage that it provides information not 
only about the shape of the correlator, but also through 
its normalization about the pion multiplicity distribution 
which is lost in the second definition. We will show that 
it can be exploited to search for multi-pion symmetriza- 
tion effects and may thus be a useful ingredient in the 
HBT analysis of 2-pion correlation functions. 



II. A SIMPLE EXAMPLE 

To illustrate the importance of the normalization of 
the 2-pion correlator let us start with a simple example 
in which we consider the following multi-pion states: 



/))„, = A,„exp[(Bt)™]|0) 
= j d^pj{p)a^ip). 



(8) 
(9) 



The states (g) are normalizable for m < 2; the norma- 
lization Am ensures m{4>\<^m = 1- For m = 1, is a 
standard coherent state ||7|] with Ai — exp(— no/2) and 
no = J d^p\j{p)\'^. Then the pion multiplicity distribu- 
tion is of Poisson form, 



Pin) = -7exp(-no), 



(10) 



and the two-pion correlators (^ and (^) are given by 

C'{pi,P2)^C"{pi,P2) = l. (11) 

For m = 2 we have A2 = (1 — 4n§)3, and the pion mul- 
tiplicity distribution is given by 



1 



n = 2fc + 1 odd, 
n = 2k even. 



Calculating the correlation functions we find 



1 

4^' 



(12) 
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One observes that now is different from C^^ due to the 
fact that the pion multiplicity distribution is no longer 
of Poisson form. However, although |0)2 is clearly not a 
coherent state, C^^ is again equal to 1. The use of the 
second definition (||) thus does not allow to distinguish 
between the states \4>i) and |02}, whereas the first defini- 
tion (^) clearly does. One may, of course, argue that the 
only difference between and C^^ is the normalization 
factor which can be obtained independently by measur- 
ing the pion multiplicity distribution. Our point is that 
important information on the pion multiplicity distribu- 
tion can also be extracted directly from the properly nor- 
malized correlation function, and that this opportunity 
should not be given away by working with probabilities 
rather than directly with the measured cross sections. 



III. MULTI-BOSON EFFECTS ON THE 
CORRELATOR AND ITS NORMALIZATION 

We will now consider a more physical model and show 
again that the use of the second definition leads to a loss 
of interesting information about the source. It is well 
known that in relativistic heavy-ion collisions the pion 
multiplicity is so large that it may be necessary to take 
multi-pion BE correlations into account P,p|,p"5|-p5[ . In 
the following we will use a specific class of density matri- 
ces for multi-pion systems p6[-p^ to study multi-pion BE 
correlation effects on the two-pion correlation function, 
thereby generalizing the conclusions of Ref. 10 . For this 
class of ensembles it was shown in [|6| Q that, after in- 
cluding multi-pion correlation effects, the two-pion and 
single-pion inclusive distributions can, in the notation of 
^^|-|l8|, be written in the following simple form ]29|] : 

N,ip)^EpHip,p), (14) 
iV2(pi,P2) = Ep,Ep, [H{pi,pi)Hip2,pi) 

+H{pi,p2)Hip2,Pi)] , (15) 

OO 

i?(pi,P2) = ^G,(pi,P2). (16) 

1=1 

The Gi (p, q) are defined as 

(17) 



where p{Pi,Pj) is a Fourier transform of the source emis- 
sion function g{x, K): 



p{Pi,Pj) = J d'^xg{x,Kij) e" 



(18) 



Here Kij={p^ +Pj)/2 'Andq^^pi -p^, q%_^Ep^ - j^p^ 



Inserting the expressions (|14|,|15|) into Eq 



Ep 

one obtains 



rl(r, r,\-^ ^ H{pi,p2)H{p2,Pi) 
H{pi,Pi)H{p2,P2) 



(19) 



This correlator goes to 1 as q cxd and to 2 as q — > 0. 
(Final state interactions are neglected here.) Thus even 
dramatic multi-boson effects as discussed below do not 
affect the intercept of the correlator — although 
they change the multiplicity distribution towards Bose- 
Einstein form they do not lead to genuine phase coher- 
ence. 

Explicit expressions for the pion multiplicity distri- 
bution and its first two moments {n), {n{n — 1)) for 
the model studied here can be found in |16|-p9[. Since 
H{pi,P2) = H*{p2,Pi), the second term in Eq. (|9|) is 
always positive, ensuring that C — 1 is positive definite. 
The normalization conditions (H) therefore imply that for 
the class of systems studied here and in jlGHTot one has 
always (n(ri — 1)) > (n)'^ (see Fig. 3 below). Obviously, 
Eq. (|l^) can therefore not apply to systems with multi- 
plicity distributions P{n) which give {n{n — 1)) < (n)^ 
(e.g. for systems with fixed event multiplicity p^ , |20| , p5[ | ) . 

The structure of ( pj| ) permits to introduce, in analogy 
to (u8), a modified source distribution S{x,K) via 



H{pi,P2) = J d^xS{x,K) e^9- 



(20) 



such that the correlator ( [T9| ) can be written in the form 
(^. S{x, K) is related to the original source distribution 
g{x,K) via Eqs. (p^)-(p^). It includes all higher order 
multiparticle BE symmetrization effects. When inter- 
preting measured single particle spectra and two-particle 
correlations one must keep in mind that the extracted 
information on the source corresponds to the effective 
source distribution S{x,p) rather than to the emission 
function g{x,p). The following example shows that these 
two functions can differ considerably; but we will also see 
that an important clue as to how much they differ will 
be provided by the normalization of the correlator. 

As shown i n |l^ p^ the recursion relations for the func- 
tions Gi in (|l7) can be solved analytically for the class 
of model ensembles studied here if the following source 
distribution g{x,p) is assumed: 



9{r,t,p) = no 
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g{r,t,p) is the Wigner density of the source in the ab- 
sence of multi-pion symmetrization effects. It is obtained 
in by folding the Wigner densities of individual 

Gaussian wavepackets with a classical phase-space dis- 
tribution jOciass for their centers (Eq. (19) in see 
also [^). The parameters R,A in ( pi] ) arc thus com- 
binations of the wavepacket width a with the spatial 
and momentum space widths i?ciass and Adass of Pciass 
(see Eqs. (20,21) in (28j). While the width parameters 
i?ciass, Aciass of the classical distribution pciass are uncon- 
strained, the widths R, A of the Wigner density g{r, t, p) 
which result from the folding procedure always satisfy 
the quantum mechanical uncertainty relation i?A > fi/2. 

The input multiplicity distribution is taken to be Pois- 
sonian as in (^0|); its mean value no can be interpreted 
as the mean pion multiplicity in the absence of Bose- 
Einstein correlations ||l^,l^. By inspection of Eqs. (|l6|)- 
( |l8| ) one easily convinces oneself that the instantaneous 
character of (|2l] ) carries over to the effective source dis- 
tribution SYx, p). Using the analytical expressions from 
Refs. [^-19 we compute 
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EpH{p,p)^EpJ d^xSix,p) 
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FIG. 1. Multi-pion correlation effects on the single-pion 
spectrum. The dash-dotted line corresponds to the input dis- 
tribution / d'^x g{x,p) with i? = 3 fm and A = 140 MeV. The 
other lines correspond to the measured distribution for vari- 
ous values of the average pion multiplicity (n) per event. Also 
given are the corresponding average phase space densities d 
(see text). 

as well as the normalized single-pion probability distri- 
bution in momentum space 

oo 



Pi''\p) = ^Jd^xS{x,p) = 
The mean pion multiplicity (n) is given by 



(n) = / d?pd'^xS{x,p) = 



i=l 



Jd'pG,{p,p). (24) 



-Pi"^ (p) is a function oi p = \p\ only. 
It is shown in Fig. 1 for different observed average pion 
multiplicities (n). Next to the value (n) we also give the 
average pion phase-space density of the system, 



in) 



(2i?A) 



3 " 



(25) 



One sees that as d increases the pions concentrate in 
momentum space at low momenta. This reflects their 
bosonic nature: pions like to be in the same state. 

The instantaneous nature of the (effective) emission 
functions g{x,p) and S{x,p) (see (pl|)) allows for inver- 
sion of the Fourier transform (|20|): writing S{x,p) = 
S{r,p) S{t) we have 



S{r,K) 



d\ 
(2^ 



H(K 



2 K 

2 ' 



f) 
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We define the normalized source distribution in coor- 
dinate space via 



Jd^KS{r,K) 
J d^Kd^rS{r, K) 
1 fd^Kd^q 
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FIG. 2. Multi-pion correlation effects on the spatial pion 
distribution. The dash-dotted line corresponds to the in- 
put momentum distribution J dtd^p g{x,p) with parameters 
i? = 3 fm and A = 140 MeV. The other lines correspond 
to the effective spatial distribution which would be extracted 
from HBT measurements, for various values of the average 
pion multiplicity (n) per event. Also given are the corre- 
sponding average phase space densities d (see text). 
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it is a function of 
7f) is known ana- 



Due to the spherical symmetry of ( 
r = |r| only. The function H in (26 
lytically ]l^-p^ to be a simple Gaussian in q, rendering 
the Fourier transform trivial. The resulting P^") (r) is 
shown in Fig. 2 for different pion phase-space densities. 
One sees that with increasing phase-space density the 
multi-pion BE correlations also lead to a concentration 



3 



of the pions in coordinate space. The fact that multi- 
pion BE correlations lead to a reduction of the HBT ra- 
dius has been observed previously ||l^,0,|5| . The radius 
extracted from HBT interferometry reflects the typical 
length scale of P^"^(r); it is always smaller than the in- 
put geometric radius R of the source and depends on the 
mean pion multiplicity per event. 



IV. THE HIGH DENSITY LIMIT 



Taking the limit of a highly condensed Bose gas, d = 
f^2iiAy^ ~* °° ' multiplicity distribution and 1- and 
2-particle spectra can be determined analytically jlj]: 



P{n) 



(n) 



N2{pi,P2) = 2 Niipi) Ni{p2 



exp 



2A2 



A 
2R 



< A^ 



In this limit the correlation functions are 

C^(P1,P2)=2, C''{pi,P2) = 



1. 



(28) 

(29) 
(30) 
(31) 

(32) 



Multi-pion BE correlations change the original Poisson 
distribution into the Bose-Einstein multiplicity distribu- 
tion (p8|). Correspondingly, ^"'"^yP^ changes from 1 to 2. 
In (|^) this change exactly compensates the fact that the 
correlator no longer decays as a function of q=pi P2, 
and from the resulting C^^ = 1 one might thus be misled 
to conclude (incorrectly) that the source exhibits phase 
coherence. 

In Fig. 3 we show the ratio i^n{n-i)) ^ function of 
the average pion phase-space density. In the lower dia- 
gram we plot it as a function of the ratio of input pa- 
rameters 'no/{2RA)^, in the upper diagram we use as a 
measure of the phase-space density the analogous ratio 
formed with the measured average multiplicity (n). For 
low phase space densities and large systems (2i?A ^ 1) 
(Tt("-i)) close to 1; BE symmetrization effects on the 
observed multiplicity distribution are then negligible, and 
it becomes equal to the Poissonian input distribution 



with {n{n - 1)) = {ny = 



For smaller systems 



(2i?A ~ 1) the observed multiplicity distribution be- 



comes non-Poissonian, with 



< 1, even in the 



{n(«-l)> 

limit of vanishing multiplicity, d = (n) / {2RA)^ 0: for 
the model discussed here one finds analytically 
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no^O 



("("-1)) 

(n}2 



= 1 



1 



(2i?A)= 



< 2. 



(33) 
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FIG. 3. Multi-pion correlation effects on the normaliza- 
tion factor as a function of the average pion phase-space den- 
sity, for different total phase space volumes 2RA — 1, 2, 4, 10 
as indicated. The corresponding input parameters were 
R = 0.7, 1.4, 3, and 7 fm, respectively, as well as A = 140 
MeV in all cases. For details see text. 

For large phase-space densities d>l, the ratio ^„(„_]^)^ 
decreases, eventually approaching for d — > oo the value 
0.5 which reflects Bose-Einstein statistics. The critical 
phase-space density for the transition from Poisson statis- 
tics with {n{n — 1)) = (n)^ to Bose-Einstein statistics 
with {n{n — l)) = 2{n)'^ depends on the total phase-space 
volume (2i?A)3, but for 2i?A > 1 it occurs near d ~ 0.3. 
For d ^ 1, multi-boson symmetrization effects become 
dominant. The Bose condensation limit is reached at a 
finite critical value for the mean input multiplicity tiq: 

\3 



no 



{RA 



> 1 



(34) 



As uq approaches the critical value (which for large sys- 
tems 2RA > 1 corresponds to nJ{2RA)^ « i), the 
observed mean multiplicity (n) and phase-space density 
d (as well as the total energy) go to infinity jl6j . In this 
sense the limit uq — s- Uc here is analogous to the limit 
/ijr — > TOtt in a thermalized pion gas of infinite volume in 
the grand canonical formalism. 
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V. NORMALIZED CORRELATION FUNCTIONS 
FROM EXPERIMENT 

A direct experimental determination of the phase- 
space density d = {n)/{2RA)^ in high energy coUisions 
is not easy. The sources created in such cohisions feature 
strong collective expansion , and therefore only small 
fractions of the total collision region (so-called regions of 
homogeneity) contribute effectively to the two-particle 
correlator ||^. This means that in the parametrization 
@ we should use i?^ = j^2^^^ + 1/(4A2) where i?hom 
is the (pair momentum dependent) homogeneity radius 
which, in the absence of strong multi-pion effects, is equal 
to the HBT radius parameter i?HBT- Without more de- 
tailed model studies it is then, however, unclear what 
fraction of the total observed multiplicity (n) comes from 
a single such homogeneity region. 

On the other hand Fig. 3 suggests that, within our 
model class of event ensembles, the ratio , ... is a 
useful indicator for the average phase-space density d in 
the source and thereby also for the expected multi-pion 
symmetrization effects on the 1- and 2-particle spectra 
which need to be taken into account in an extraction of 
the source size from HBT measurements. 

In the experiment one usually fits the two-particle cor- 
relator with the functional form 

C^^P(P1,P2) = C<=''P(q,K) = A/-(l+A|/(q,K)p) , (35) 

where the function / vanishes as q — > oo. Obviously, 
TV depends on the chosen definition of the correlation 
function: For the definition (^) the normalization is al- 
ways TV^ — 1 (see Eq. dl^)), while for the definition ^ 

it is TV^^ = (n("^i)) (^^"^ ®)- ^^^^ both cases 
Af is well-defined and thus should not be treated as a 
free fit parameter. Therefore we now discuss shortly an 
algorithm for the experimental construction of the two- 
particle correlator which is guaranteed ||3^ to give the 
correct value for Af, without relying on an actual mea- 
surement of the multiplicity distribution. 

We write the single-pion inclusive distribution as 



Er. 



We- 



(36) 



These definitions satisfy the normalization conditions (||). 
While N2{pi,P2) is constructed by selecting pion pairs 
from the same events, the denominator iVi(pi)A^i(p2) 
can be generated by combining pion pairs from different 
events ]33|,p|- p^p^ -p^] . The proper prescription is 



Niipi)Niip2) 



E E 

iVev(iVev - 1) 



JVov 

E 

i.j=l 



Vt,3{Pl,P2) (38) 



where P'i.j{pi,P2) = i'i{Pi)i'j{P2)- One easily checks that 
d^pi d^p2 



E E 



■Ni{pi)Ni{p2) = {nY 



(39) 



The ratio of (^7|) and (|3^) thus gives the properly nor- 
malized correlator . The above equations are true for 
unbiased events. Trigger biases and limited experimen- 
tal acceptances can induce residual correlations in the 
event-mixed "background" (^sj) which must be corrected 
for separately (see |^ for an extensive discussion). 

For large N^v the evaluation of (^) is very time con- 
suming; it also leads to a statistically unnecessarily accu- 
rate result for the denominator in (^. In practice one can 
live with fewer event pairs for event mixing, by replac- 
ing in (|3^) the number iVov by a much smaller number 
N'^^. As long as Nl,^{N'^^ — 1) TVov one can still ensure 
that the contribution of the denominator to the statis- 
tical error of the final correlation function is negligible 

H- 

The correlator C^^ differs from only by the different 
normalization. It can be constructed by taking the ratio 
of the following two expressions B : 



Pl{Pl)Pl{P2) 



E E 



P2{pi,P2) = E ^>m(Pi,P2) , (40) 

pairs 



TV" , 

pairs i,i=l 



51 i^M(Pi.P2). (41) 



N';^^^^ and iV" • „ are the total numbers of "correlated" 
and "uncorrelated" pion pairs, respectively: 

A^pairs = ^cv • {n{n ~ 1)) , (42) 



4=1 



^pairs = ^cv(TVcv - 1) ' {nf 



(43) 



iVcv is the total number of collision events, and I'iip) is 
the number of pions with momentum p in collision i. The 
two-particle distribution can be expressed as 



iV2(Pl,P2) = %^E'^M(Pl,P2). 



(37) 



where Vi,i{pi,P2) is the number of pion pairs with mo- 
menta (pi,P2) in collision event i, and the double index 
i indicates that both particles are from the same event. 



VI. CONCLUSIONS 

We have shown that in principle the normalization of 
the two-particle Bose-Einstein correlation function con- 
tains valuable information on the the multiplicity dis- 
tribution of the event ensemble. Both theoretically and 
experimentally the absolute normalization of the corre- 
lation function should thus be controlled as well as pos- 
sible. We presented a variant of a previously suggested 
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experimental algorithm [Q for the construction of the 
correlator which guarantees correctly normalized corre- 
lators. Within a specific model class of event ensembles 
which recently received extensive theoretical attention we 
showed that in systems with large pion phase-space den- 
sities multi-pion symmetrization effects can lead to inter- 
esting measurable effects on the normalization of the cor- 
relator. We suggest a careful study of this normalization 
as an alternate method for searching for strong multi- 
pion symmetrization effects in high-multiplicity hadronic 
and heavy-ion collisions. 
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